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We present here the system size dependence of balance energy for semi-central and pe- 
ripheral collisions using quantum molecular dynamics model. For this study, the reactions 
c5 ■ of Ne^"+Ne^°, Ca^o+Ca^o, Ni^'^+Ni^'^, Nb'^^+Nb^^, Xe"i+Xei3i and Au^'^^+Au^^^ are 



simulated at different incident energies and impact parameters. A hard equation of state 
along with nucleon-nucleon cross-sections between 40 - 55 mb explains the data nicely. In- 
terestingly, balance energy follows a power law oc A'^ for the mass dependence at all colliding 
geometries. The power factor r is close to — ^ in central collisions whereas it is — | for 
^ ^ peripheral collisions suggesting stronger system size dependence at peripheral geometries. 

^ , This also suggests that in the absence of momentum dependent interactions. Coulomb's in- 



teraction plays an exceedingly significant role. These results are further analyzed for nuclear 
dynamics at the balance point. 
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I. INTRODUCTION 



I Heavy-ion collision at intermediate energies is a complex phenomenon that depends crucially on 

. the interplay of the mean field and nucleon-nucleon binary collisions. The dominance of either of 

these ingredients can shed light on the properties of nuclear matter at the extreme of temperature 
and density. Among various observables and non-observables in heavy- ion collisions, directed 
collective flow is considered to be a very important observable due to its extreme sensitivity towards 
the model ingredients and physics involved 1-^. The very same observable has also been robust 
for the understanding of nuclear equation of state - a question that is still open and extensively 
debated in the literature j^-lO]. The directed collective flow has been reported by many authors 



to be positive at higher incident energies whereas same is deeply negative at low incident energies. 
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While going from the low to high incident energies, there comes a typical point where no preference 
is observed for the flow, therefore, it disappears. This particular incident energy where attractive 
and repulsive interactions balance each other (and flow disappears) is termed as balance energy 
or Ebai [9, 
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181 ]. At the balance point, scatterings due to the attractive mean field occurring at 



negative angles are counter balancec 
occur at positive angles 
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151, 



191, 



by the repulsive action of the binary nucleon-nucleon collisions 



2C|. 



The composite dependence of the Ei,ai on the mean field equation of state (eos) and nucleon- 
nucleon cross-section (ex at at) can be sorted out by noticing the sensitivity of Ef^ai on the system 
size as well as on the impact parameter of the reaction. In this connection, very recently, Puri 
and Sood 3], conducted a very detailed and exhaustive study on the energy of vanishing flow 
over entire periodic table with masses between 20 and 394. Their study had shed light on various 
aspects of nuclear dynamics. Unfortunately, this study along with all other studies reported in 
the literature are limited to central/semi-central collisions only [19]. A very little attention has 
seen paid in the literature for semi-central and/or peripheral collisions, where physics is derived by 



the low density dynamics 



IC 



I9I423I]. Among various attempts involving non-central impact 



parameters, Magestro et. al. [17'] reported the results of disappearance of flow over large impact 
parameters. Their study, however, was limited to heavier systems only. One should keep in 
the mind that the dynamics in lighter nuclei is entirely different than that of heavier colliding 
nuclei. It remains, therefore, a challenging question to see how directed transverse flow behaves at 
perpheral/non-central collisions. We plan to address this question in this study. 



In an another study, Puri et.al. 



131 ] extracted the strength of nucleon-nucleon cross-section for 



the reaction of Zn^'^ + At^"^ at different collision geometries. They firmly indicated that for central 
collisions, a cross-section of 40 mb was good enough, whereas an enhanced value of cross-section 
was needed for peripheral collisions. It still remains to be seen whether this observation holds good 
over entire periodic table or not. 

Here we plan to understand the mass dependence of disappearance of flow in these low excited 
geometries. We shall also attempt to parameterize the balance energy for the system size effects. 
The present study is conducted within the framework of quantum molecular dynamics (QMD) 
model, which is discussed in section II. In section III, we discuss the results and summary is 
presented in section IV. 
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II. THE MODEL 



In QMD model 24l426l|. each nucleon, represented by a Gaussian distribution, propagates with 



classical equations of motion as: 



where Hamiltonian H is given by : 



dt ~ dpi' ^ ' 



dpi ^ _dH^ 

dt dvi ' ^ ' 



i 



with 



'Y'tot -yloc _|_ -y^Yuk _|_ -y^Coul 

Here V^°'^, V^^^ and represent, respectively, the local two- and three-body Skyrme, Yukawa 

and Coulomb interactions. Under the local density approximation, Skyrme part of the nucleon 
interaction can be written as: 

y'- = a/2j;p. + -A_J^,J. (5) 

i i 

Here pi is the nucleon density. The coefficients a and /3 give proper ground state properties and 
7 is used to generate different compressibilities. The different values of compressibility (7) give 
possibility to look for the role of different equations of state termed as soft (7 = 200 MeV) and stiff 
(7 = 380 MeV) equations of state. The values of a and /3 for soft and hard equation of state are 
-0.356 GeV, 0.303 GeV and -0.124 GeV, 0.0705 GeV [24]. Following many studies [l3,Q,Q,[3 
listed in the literature, we also use a hard equation of state. 

When propagating nucleons come too close to each other, they can collide elastically or inelas- 
tically depending upon the available center of mass energy. The collision probability depends on 
the cross section ctnn- Several different forms of a^N are available in the literature. In many 
studies Q, [3], one has taken constant value of aNN- Generally it is assumed to be between 20 mb 
and 55 mb. The strength of nucleon- nucleon cross-section will be extracted by simulating various 
reactions with between 40 and 55 mb. The above QMD model has been shown to be very 



useful in many observables /non-observables [25|, [26 1 
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FIG. 1: Average < -pxl^ > (MeV/c) as a function of scaled rapidity ycm/Ubeam- We display the results 
(anticlockwise) at impact parameters b = 1, 2, 3, 4, 5, 6 and 7 fm. 

III. RESULTS AND DISCUSSION 

For the present study, symmetric reactions of A^^e^" + ATe^", Ca^° + Ca^°, Ni^^ + Ni^^, Nb^^ + 

are taken. The entire colhding geometry is considered with 
main emphasis on the scmi-ccntral as well as on the peripheral collisions. Throughout the present 
analysis, a hard equation of state along with ajviv = 40- 55 mb has been used. In fig.l, we display 
the changes in the transverse momentum < Px/A > as a function of the rapidity distribution 
ycm/ybeam- Here we vary impact parameter from b = 1 to b = 7 fm at a step of unity. We see 
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FIG. 2: Same as in fig. 1 but here we display the results (anticlockwise) at energies 40, 60, 80, 100, 150, 
200 and 400 MeV/nucleon. 

that for b = 1 fm, flow is nearly isotropic and slope is zero. The negative slope increases as we 
move from central to peripheral collisions. This demonstrates that for larger impact parameters, 
matter is attractive in the presence of reduced number of collisions. It is worth that the slope of 
the < Px/-^ > depends crucially on the incident energy. In fig. 2, we display the < Px/A > for the 
same reaction as fig. 1 but at b = 4 fm and for different incident energies. We see that slope is 
negative at small incident energies whereas it turns less negative with the increase in the incident 
energies. At particular energy (between 100 and 150 MeV/nucleon), we see a change in the slope 



6 




50 100 150 200 250 300 50 100 150 200 250 300 



Time (fm/c) 

FIG. 3: The time evolution of total directed transverse momentum < P^*'' > along with mean field and 
collision parts at the energy of vanishing fiow. Solid symbols represent Ca"^^ whereas open symbols are for 



from negative to positive values. This slope is steeper at higher incident energies. 

In Fig. 3, the transverse directed momentum, < P^*** > is plotted as a function of reaction 
time for the collisions of Ca'^^+Ca^^ and Au^^"^ +Au^^''^ using a^N = 40 and 55 mb, respectively. 
The displayed results are at reduced impact parameters b/bmax (b) = 0.45 and 0.60. The plots 
are displayed at the energy of vanishing flow i.e. balance energy. For a better understanding, we 
divide the total transverse momentum into contribution resulting from the mean field and from 
the collision parts. From the figure, one observes that the transverse momentum due to mean field 
(representing squares) is being equated by that due to binary nucleon-nucleon collisions (circles), 
resulting in net zero fiow (triangles) at the end of the reaction. The mean field contribution 
dominates the reaction during early evolution resulting in negative (attractive) transverse fiow. 
It is worth mentioning that if mean field dominates, the transverse flow will be negative whereas 
the dominance of nucleon-nucleon collisions will result in to a positive fiow. At ultra low incident 
energies, the phenomena emerging due to mean field alone are dominant j^]. The reaction persists 
longer for heavier systems (up to 100 fm/c or so for Au^^"^) compared to lighter colliding nuclei 
(which is about 50 fm/c). This happens due to lower energy of vanishing of fiow in heavier systems 
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FIG. 4: The total binary collisions (< A^g^^ >) at 200 fm/c as a function of reduced impact parameter (b). 
Lines are drawn to guide the eye. Here the reactions of Ne'^^+Ne'^^, Ca^^+Ca'^^, Ni^^+Ni^^, Xe^^'^+Xe^^'^ 
and Au^^"^ +Au^^'^ are taken. 

compared to light systems. The different values of the flow due to mean field using different a^N 
happens due to different energies of vanishing flow. 

We plot in Fig.4, the number of allowed collisions < Ng°}/^ >, for various systems as a function 
of reduced impact parameter b. This dependence of binary collisions on impact parameter yields 
several interesting aspects: For the lighter masses, there is a marked enhancement in the binary 
collisions at peripheral collisions. This happens due to exceedingly high balance energy at these 
impact parameters in lighter colliding nuclei. For example, the E^ai for Ne^'^+Ne^^ at 6 = 0.15 
is 76 MeV/nucleon, which goes up to 360 MeV/nucleon for h = 0.60. On the other hand, binary 
collisions for heavier colliding nuclei (like Av}^"^ + Av}^"^) are almost independent of the impact 
parameter. This is due to the fact that the energy of vanishing flow is very small in Av}^"^ +Av}^'^ 
reaction (between 40 and 78 MeV/nucleon) and at these low incident energies, most of the binary 
collisions are blocked, and as a result no significant impact parameter dependence is obtained for 
binary collisions in heavy-ion systems. 

In figs. 5 and 6, we plot the energy of vanishing flow as a function of total system mass A (mass 
of projectile + mass of target) for entire mass range between 40 and 394 units. Remarkably, our 
present calculations with unn between 40 and 55 mb explain experimental data for h = 0.45 nicely, 
whereas some deviation appears in the case of 6 = 0.60. The data shown in the figure are for the 
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FIG. 5: The balance energy as a function of system mass at b — 0.45. The data points are taken from 
references mentioned in text. Shaded region represents the energy of vanishing flow obtained with aNN = 
40 and 55 mb 



reactions of ^r^o + Sc^^Ca'^^ + Ca^o, Nj^^ ±_Nj^ ,Fe^^ 
TVi^s, Zn64 + Ti^s and Au^^^ + Au^^'^ 



171, 



2QI, 



21 



23 



^g58f^^86 ^ Nb^^,Ar^^ + Af^, Zn^^ + 



291 ]. The scarcity of data at larger 



impact parameters restricts the similar detailed analysis. We notice that different cross-sections 
have larger influence for light colliding nuclei that reduces significantly for heavier systems. The 
cause suggests the fact that the balance energy for heavier system is much smaller compared to 
lighter nuclei. As a result, cross-section plays an insignificant role in heavier systems. This mass 
dependence can be parameterized by a power law of the form cx A'^ . As we see, the r for b = 
0.45 is -0.44 and -0.37 for our calculations with aNN = 40 and 55 mb respectively, whereas it is 
-0.66 and -0.57 for = 40 and 55 mb at 6 = 0.60. The deviation of our calculated energy 

of vanishing flow from experimental data suggests that the dynamics of low energy domain needs 
further investigation. We also notice that around A = 100 units, no clear mass dependence appears. 
This could also be due to the fact that impact parameter in these measurements is not fixed. It 
varies with reaction. The impact parameter variation can also have effect on the energy of vanishing 
flow. 

Fig. 7 shows our calculated balance energy versus system mass at reduced impact parameters 
of 0.15, 0.30, 0.45 and 0.60. All calculated points in the energy scale are fitted with a power law 
oc A"^ . One notice that the slope of plots increases with impact parameter. The values of power 
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FIG. 6: Same as fig.5, but at S = 0.60. 




FIG. 7: The balance energy Ei,ai against system mass at different reduced impact parameters. The displayed 
results are using QMD model with hard EOS at 55 mb cross-section. 

factor T are -0.25, -0.22, -0.34 and -0.53 respectively, for b = 0.15, 0.30, 0.45 and 0.60 (Note that 
in this figure and fitting, Ne^^+Ne^^ reaction has been excluded, therefore slope deviates shghtly 
compared to fig.6).We see that the value of r depends strongly on the colliding geometry. For 
instance, for central collisions, E^^ai is 50-57 MeV/nucleon for Ca^° -|- Ca^ reaction, whereas it 
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FIG. 8: The time zone for which density remains higher than half of the normal matter versus system size 
for different impact parameters. 

increases to 149.5 MeV/nucleon for b = 0.60. 

For central collisions, we observe that the power factor is close to (-1/3). This is a well known 



trend reported by number of authors earlier 
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Q, 13] and has been very well justified in the 



literature in terms of interplay between the mean field and binary collisions. The increase in the 
slope of Eijai with impact parameter can be due to failure of static equation of state to reproduce 
transverse motion of particles. The effect is more severe for lighter nuclei. This is partially due to 
the fact that with weaker binary nucleon-nucleon collisions, Coulomb's repulsion plays a decisive 
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FIG. 9: Maximal time for collision rate and average density as function of system size. 

role. This role is very strong in heavier nuclei and at peripheral collisions, therefore making slope 
steeper. It still remains to be seen how momentum dependent interactions can alter these findings. 

To further understand this similarity in the balance energy obtained at different impact pa- 
rameters, we display in fig.8, the evolution time for which density remains higher than half of the 
matter density i.e. the time for which p > O.Spo- This gives an important information about the 
interaction time at different impact parameters. We display the results at 6 = 0.15, b = 0.3, b = 
0.45, b = 0.6, b = 0.75 and b = 0.9. Interestingly, we see that in all cases, time zone for different 
masses can be explained with the help of a power law oc c'A^' . The factor t' increases systemati- 
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cally from central to peripheral reactions. It is 0.42 for central whereas it increases to 0.67 for b = 
0.9. A careful look reveals that the effect is less for heavier nuclei (for Au+Au it is 80 fm/c for b 
= 0.15 whereas, it is 20 at 6 = 0.9). On the other hand, it has drastic effect for lighter one where 
it slashes to 10% of central values. This is the cause why balance energy changes drastically for 

lighter nuclei compared to heavier ones and we have different slopes compared to peripheral ones. 

In fig. 9, we display the maximal time of average density and collision rate for the same systems 
as shown in fig.8. Here we also see again a similar trend as seen in fig.8. There is a monotonous 
increase in the slope from central to peripheral collisions. 

IV. SUMMARY 

In the present work, we analyzed the disappearance of flow (or alternately the balance energy) 
for semi-central and peripheral collisions using quantum molecular dynamics model. The present 
study was undertaken by using a hard equation of state along with constant cross-sections between 
40 and 55 mb. We sec a mass dependence oc A^^^^ for central collisions, whereas this dependence 
becomes steeper with impact parameter. This deviation from -1/3 trend suggests exceeding role 
of Coulomb's repulsion in low density and excited situations. Nuclear dynamics at the balance 
point indicates drastic changes in lighter system compared to heavy one in agreement with balance 
energy. 

This work is supported by Department of Science and Technology, Government of India. 
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